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We  have  performed  molecular  dynamics  calculations  on  ClXHjOn  and  FCHgOln 
(n=2,3,...15)  clusters.  The  calculations  show  that  the  F  ion  is  solvated  in  these  clusters,  while 
Cr  remains  attached  to  the  water  in  the  clusters.  We  also  obtained  the  minimum  energy 
structm-es  for  the  Cr(H20)„  and  F'dijOln  (n= 6,7,8)  clusters.  From  the  comparison  of  these 
structures  with  the  dynamical  structvu-es  we  conclude  that  the  solvation  of  F'  ion  is  due  to  the 
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Abstract: 

We  have  performed  molecular  dynamics  calculations  on  ClXHgOln  and  F  (H20)n 
(n=2,3,...15)  clusters.  The  calculations  show  that  the  F'  ion  is  solvated  in  these  clusters,  while 
Cr  remains  attached  to  the  water  in  the  clusters.  We  also  obtained  the  minimum  energy 
structures  for  the  ClXHjOln  and  FXH20)„  (n= 6,7,8)  clusters.  From  the  comparison  of  these 
structures  with  the  dynamical  structures  we  conclude  that  the  solvation  of  F  ion  is  due  to  the 
entropy  effect. 
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I.  INTRODUCTION 

The  study  of  hydration  of  ions  by  few  water  molecules  provides  important  information 
for  oxir  understanding  of  solvation  phenomena.  This  study  received  a  new  impetus  due  to 
recent  experimental  and  theoretical  work  on  this  subject  [1-23].  The  first  experimental  data 
on  these  systems  were  obtained  by  Kebarle  and  his  coworkers  [1,2].  Using  mass  spectroscopic 
techniques  they  measured  the  temperatme  dependence  of  the  equilibrium  constants  K  for 
the  reaction: 


Once  the  temperature  dependence  of  the  equilibrium  constant  is  known,  the  heat  of  the 
reaction  (1)  can  be  obtained  from  the  Van’t  Hoff  equation: 


HH 


0 


3(1/7) 


(2) 


In  its  turn  the  knowledge  of  the  heat  of  the  reaction  (1)  can  be  used  to  test  the  quality  of  the 
interaction  potentials  between  the  ion  and  water  molecules.  Unfortunately  the  experimental 
measurements  of  the  heats  of  the  r  eaction  (1)  could  be  performed  only  for  very  small  sized 
clusters  (with  up  to  n  =  8  water  molecules  ),  thus  limiting  our  testing  abilities  of  the  quality 
of  the  potential.  Moreover,  from  the  knowledge  of  the  heat  of  the  reaction  it  is  hard  to 
extract  information  on  the  structure  of  the  clusters.  Recently  an  improvement  in  the 
experimental  situation  has  occurred  due  to  the  application  of  photoelectron  spectroscopy 
techniques  to  anion-water  clusters  [10].  Prior  to  its  use  on  the  halide-water  clusters 
photoelectron  spectroscopy  demonstrated  its  value  in  the  study  of  the  electron  solvation  in 
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water  clusters  where  it  allowed  one  to  perform  a  comparison  between  calculations  and 
experiment[24,25]. 

As  of  today  we  have  available  to  us  data  from  the  photodetachment  spectra  of  the 
following  clusters:  Cr(H20)„  (n=  BrXHjOln  (n=  1,...15)  and  r(H20)„  (n=  [26].  Very 

recently  we  have  demonstrated  that  the  experimental  stabilization  energies  for  I  ,  and 
Br  and  a  calculated  stabilization  energy  for  Cl'  in  these  clusters  satisfy  the  relationship 

-  CAn)  (3) 

where  f(n)  is  a  universal  function  that  depends  only  on  the  number  of  water  molecules  in  the 
cluster  and  C  is  the  proportionality  constant  [27].  The  existence  of  the  relationship  described 
by  equation  (3)  means  that  the  structures  of  calculated  CrdijO)^  and  experimental  X  (H20)„ 
(n=2,.,.15,  X=Br,  I )  clusters  are  similar  [28].  And  since  we  observed  in  our  simulations  that 
Cl'  ion  is  located  on  the  surface  of  these  water  clusters,  from  the  similarity  argument  we 
therefore  concluded  that  Br'  and  I'  ions  are  not  solvated  by  water  clusters  either. 

As  far  as  the  clusters  with  Cl'  are  concerned  we  observe  the  following: 

a)  Our  calculated  stabilization  energies  for  C1'(H20)„  (n  =  2...7)  clusters  are  in  a  good 
agreement  with  the  experiment. 

b)  The  experimental  £,^3^  (for  all  the  available  experimental  data)  in  Cl'(H20)n  (n=l,2,...7) 
clusters  fit  the  universal  curve. 

c)  The  quantum  mechanical  calculations  performed  on  Cl(H20)i4  cluster  indicate  the 
preference  for  the  Cl'  ion  to  be  on  the  surface  of  the  water  cluster  [29]. 

All  these  factors  in  our  opinion  provide  support  to  the  previously  proposed  conjecture 
that  the  Cl'  ion  is  also  located  on  the  surface  of  the  water  clusters  which  contain  up  to  20 
water  molecules  [30,31].  Nevertheless,  the  final  judgment  on  the  location  of  Cl'  in  these 
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clusters  should  be  given  only  after  the  photodetachment  experiment  will  be  performed. 

Considering  the  previously  given  information  about  the  structure  of  hahde  amon/water 
clusters  it  is  therefore  interesting  to  find  out  which  properties  of  the  ion  influence  the 
structxure  of  these  clusters.  In  our  previous  study  of  Cl  (H20)2o  cluster  we  observed  that  the 
polarizabiUty  of  the  ion  does  not  substantially  change  the  structure  of  the  cluster,  but  the 
sign  of  the  ionic  charge  has  a  profound  influence  on  the  structure  [31].  Another  ionic 
property  that  may  have  a  dramatic  effect  on  the  cluster  £>cructure  is  the  ion  size.  To  study 
the  dependence  of  the  cl\ister  structure  on  the  ionic  size  we  could  change  the  ion-water 
potential  in  such  a  way  that  only  the  size  of  the  ion  would  change,  and  study  the  effect  of  this 
change.  Instead  we  decided  to  study  a  realistic  system  of  the  F  (H20)n  clusters  and  compare 
the  results  with  the  results  from  our  calculations  on  Cl  (H20)n  clusters.  In  this  way  we  can 
also  complete  our  study  of  structural  properties  of  small  water  clusters  (ni  15)  with  embedded 
halogen  ions. 

II  METHODS 

Potentials. 

To  describe  the  hahde  ion-water  and  water-water  interactions  we  used  the  POLl 
potential  model  which  includes  many-body  effects  [21,22,23].  In  this  model  the  water 
molecule  is  assumed  to  have  three  charges  placed  at  the  positions  of  two  hydrogens  and  the 
oxygen.  The  Lennard-Jones  center  is  placed  on  the  oxygen  site.  The  OH  distance  in  the 
model  is  1  A  and  HOH  angle  is  109.47  ° .  Three  polarizable  centers  are  placed  at  each  atomic 
site  in  water.  The  ion  in  this  model  is  considered  to  be  a  point  charge  placed  at  a  Lennard- 
Jones  center.  Additionally  a  polarization  center  is  also  placed  at  the  same  location. 

The  total  energy  of  the  system  is  given  by  the  following  expression: 
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(4) 


where 


(5) 

i<i 


is  the  sum  of  all  Coulomb  interactions  between  all  i  and  j  charge  sites  and 

-  E  (6) 

•<f 


is  the  sum  of  the  Lennard-Jones  interactions  between  a  and  0  Lennard-Jones  sites.  The 
polarization  energy  is  obtained  from  the  expression 
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where  aj  is  the  polarizability,  /t,  is  the  induced  dipole  at  the  polarizable  i-th  center  and  Tjj  is 
the  distance  vector  from  site  i  to  site  j.  The  values  of  the  induced  dipoles  at  the  polarizable 
centers  are  calculated  self-consistently  using  relationships  (8)  and  (9) 
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where  E;  is  the  total  electrostatic  field  at  the  i-th  center  and  E  °  is  the  field  due  to  permanent 
charges.  The  dipole  tensor  Ty  in  eq  (9)  is 


The  potential  energy  due  to  the  three  body  exchange  repulsion  of  the  ion-water  trimer  is 
expressed  in  terms  of  ion-oxygen  distances  r^  and  rji^,  and  the  oxygen-o:tygen  distance 

^3-bc<fy  -  5exp(-(irpexp(-Pr^)exp(-Yrp  (ID 


where  B,  /3  and  y  are  parameters.  The  potential  parameters  of  the  ’’  jLI  model  used  in  the 
present  work  are  given  in  Table  I. 

Molecular  Dynamics. 

We  used  molecular  dynamics  computer  simulations  to  investigate  the  structure  and 
dynamics  of  small  clusters  made  up  of  one  ion  and  n  water  molecules.  Equations  of  motion 
were  solved  using  the  Verlet  algorithm  [32].  The  geometries  of  the  water  molecules  at  their 
equilibrium  values  were  constrained  through  the  Shake  procedure  [33].  To  solve  self- 
consistent  equations  (8)  and  (9)  an  iterative  approach  was  used.  The  iterations  were 
continued  until  the  root  mean  square  of  the  difference  in  the  induced  dipoles  between 
successive  iterations  was  less  than  0.001  D/atom.  The  self-consistency  was  usually  achieved 
in  seven  steps  of  iteration. 

Each  cluster  was  initially  equilibrated  for  more  than  20  ps  followed  by  1  ns  of  the 
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production  rvm.  The  time  step  of  the  trjyectory  has  been  kept  at  1  fs  and  the  total  linear  and 
angular  momenta  of  the  system  were  removed.  During  the  tr^ectory  calculations  the  energy 
was  conserved  within  5-th  significant  digit. 

All  the  trfyectory  calculations  were  done  with  the  occasional  rescaling  of  the  velocities 
so  that  the  temperature  of  the  clusters  was  kept  constant.  For  the  very  small  clusters  (ni4) 
the  temperatmre  was  kept  around  275  K,  for  the  clusters  with  5:ini9  the  temperature  was 
around  250  K,  For  the  clusters  with  9-<nil5  the  average  temperature  was  dropped  to  around 
225  K.  The  temperature  of  the  clusters  was  chosen  to  be  high  enough  to  eliminate  possible 
trapping  of  the  cluster  in  a  certain  configuration  for  a  very  long  time.  We  attempted  to 
perform  some  molecular  dynamics  simulations  of  our  clusters  at  300  K,  but  observed  the 
evaporation  of  water  molecules  at  this  temperature. 

RESULTS. 

To  test  the  quality  of  the  potential  used  in  the  simulation  we  calculated  the  heats  of 
cluster  formation  for  XXHjO)^  (x=Cl,  F;  n=2,3,.--15)  clusters  and  compared  these  with  the 
experimental  data  when  these  were  available.  The  results  of  such  comparison  are  given  in 
Table  II.  As  we  can  see  from  the  Table  the  agreement  is  rather  good,  considering  the 
uncertainties  in  the  experiment  and  in  the  simulations.  To  obtain  information  on  the 
structures  of  the  clusters  we  calculated: 

a)  distributions  of  the  distances  between  the  .on  and  the  center  of  mass  of  the  water 
molecules  in  the  cluster. 

b)  ion-water  oxygen  radial  distribution  functions  (rdD  and  running  coordination  number  (cn) 
for  an  ion  in  water  clusters. 

c)  angular  distribution  functions  (adf)  for  an  ion  in  water  clusters. 

We  have  also  examined  the  snap-shots  of  cluster  configurations  obtained  from  our  runs 
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All  the  data  from  our  calculations  lead  to  the  following  conclusion:  while  the  Cl'  ion  is  not 
solvated  in  water  clusters  with  up  to  15  water  molecules,  the  F'  ion  is  solvated  in  these  water 
clusters. 

To  justify  our  conclusion,  let  us  first  consider  the  distribution  of  the  distances  from 
the  ion  to  the  center  of  mass  of  the  cluster.  We  represent  such  a  distribution  for  clusters 
with  7  water  molecules  in  Figure  la.  While  for  the  F"  ion  the  average  distance  between  the 
ion  and  the  center  of  mass  of  the  cluster  is  1.14  A,  for  Cl'  ion  this  distance  is  1.78  A.  Figure 
lb  presents  the  distribution  of  the  distances  from  the  ions  to  the  centers  of  mass  of  clusters 
with  15  water  molecules.  In  this  case  for  the  F'  ion  the  average  distance  between  the  ion  and 
the  center  of  mass  of  the  cluster  is  2.22  A,  while  for  the  Cl  ion  this  distance  is  now  3.19  A. 
Figures  la  and  lb  indicate  that  I  '  ion  is  more  solvated  by  water  clusters  than  Cl'  ion.  In 
Figure  2a  we  represent  the  ion-oxygen  radial  distribution  functions  for  clusters  with  the  F 
ion  and  in  Fig.2b  ion-oxygen  rdf  for  the  clusters  with  the  Cl'  ion.  For  the  clarity  of  the  figures 
we  consider  the  rdfs  for  the  clusters  with  3,  7,  11  and  15  water  molecules  only.  As  we  can 
see  from  the  figures,  although  some  similarity  exists  between  the  F'-oxygen  and  Cl'-oxygen 
radial  distribution  functions,  the  rdf  for  Cl'-oxygen  displays  that  a  more  pronounced  third 
solvation  shell  is  formed  around  this  ion.  As  to  the  position  of  the  first  peak  of  rdf  we 
observed  that  it  is  shifting  towards  the  larger  distances  when  the  rdf  of  F'  ion  is  considered, 
while  such  shift  becomes  insignificant  for  clusters  with  more  than  five  waters  when  Cl'  ion 
is  in  the  cluster.  This  indicates  that  the  continuous  build-up  of  the  solvation  shell  occurs 
arovmd  F',  while  around  Cl'  the  first  neighbors  shell  is  completed  rather  early  in  the  growth 
of  the  cluster.  This  also  indicates  that  the  F'  ion  is  solvated,  while  Cl'  is  not.  Additional 
confirmation  of  this  conclusion  is  provided  from  the  plots  of  running  coordination  numbers, 
depicted  in  Fig.  3a  (for  F')  and  Fig.  3b  (for  Cl').  Note  that  while  the  coordination  number  for 
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the  first  solvation  shell  is  constantly  increasing  for  F‘,  until  it  reaches  the  value  of  5.7  in  the 
cluster  with  15  waters,  the  coordination  number  for  Cl'  approaches  the  value  of  4.4  when  10 
water  molecules  are  in  the  cluster.  Direct  information  about  the  as3mimetry  of  clusters  is 
given  by  angular  distribution  functions  [31]  which  are  presented  in  Fig.  4a  and  Fig.  4b.  For 
the  cluster  where  the  ion  is  surrounded  by  waters  from  all  the  sides  the  distribution  should 
look  flat.  A  distinct  as5anmetry  in  distribution  is  observed  for  the  clusters  with  the  Cl'  ion, 
while  the  distribution  is  approaching  a  flat  character  for  the  F‘(H20)i5  cluster.  (Again,  for  the 
matter  of  clarity,  we  present  the  distribution  fvmctions  for  clusters  with  3,  7,  11  and  15  waters 
only).  Finally,  snap-shots  depicting  configurations  of  the  clusters  with  seven  water  molecules 
and  the  corresponding  ions  are  presented  in  Fig.  5a  and  Fig.  5b. 

Now,  since  we  know  that  F'  is  solvated  in  water  clusters  and  Cl'  is  not,  we  can  ask  why 
this  is  happening.  The  energy  component  analysis  which  is  given  in  Tables  Hla  and  Illb  may 
help  us  in  providing  some  clues.  For  the  energy  component  analysis  we  divided  the  total 
energy  into  the  corresponding  contributions  from  the  pair  interaction  energy  Ep^j,., 
polarization  energy  Ep^,  and  from  the  3-body  term,  so  that 

E  ~  E  +E  ,+£,  w 

tot  pair  pel  ^-body 


In  its  turn  the.  pair  interaction  energy  can  be  written  as  a  sum  of  water-water  and  water-ion 
interaction  energies;  these  can  be  divided  into  Lennar d- Jones  and  electrostatic  parts; 

E  ■  -  E  ,+E  ,,+E  .  ,+E 

pair  ww,el  wwJJ  H»,ri  wtJJ 


(13) 
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where  E**  ,,  is  the  electrostatic  part  of  the  water-water  pairwise  interaction,  u  is  the 
Lennard-Jones  part  of  the  water-water  pairwise  interaction,  E,,,-  g,  is  the  electrostatic  part  of 
the  ion-water  pairwise  interaction  and  i®  Lennard-Johes  part  of  the  water-ion 

pairwise  interaction.  As  Tables  III  show,  a  nuyor  difference  in  the  energies  for  F'  and  Cl  is 
observed  for  the  E,„j ,,  entry.  This  energy  is  lower  for  clusters  with  F'  ion,  which  may  be  due 
to  the  fact  that  the  F'  is  more  solvated.  It  may  be  also  due  to  the  fact  that  F'  ion  is  smaller 
in  size  than  Cl'  ion.  That  the  F'  is  more  solvated  than  Cl'  ion  can  be  seen  by  the  inspection 
of  another  component  of  the  energy,  E,„,  which  is  lower  for  the  clusters  with  Cl'  ion.  This 
indicates  that  the  hydrogen  bonded  network  of  water  is  less  disturbed  by  the  presence  of  Cl' 
ion,  i.e.  the  Cl  ion  is  less  solvated. 

Another  very  interesting  fact  that  can  be  observed  from  Tables  III  is  that  the  relative 
contribution  of  the  polarization  energy  into  the  total  energy  of  the  cluster  is  continuously 
diminishing  (from  -33%  for  n  =  2  to  ~  18%  for  n  =  15)  for  the  clusters  with  F'  ion,  while  for  the 
clusters  with  Cl'  ion  the  relative  contribution  of  the  polarization  energy  remains  around  the 
same  ( ~20%).  That  may  mean  that  it  is  important  to  include  the  polarization  effect  when  the 
structures  of  Cl'/water  clusters  are  calculated.  To  see  the  effect  of  the  explicit  inclusion  of 
the  polarization  into  the  calculational  scheme  we  represent  in  Figure  6  the  values  of  the 
dipole  moments  for  water  molecules  as  the  function  of  their  location  in  the  solvation  shells. 
As  we  can  see  the  value  of  the  dipole  moment  of  the  water  molecule  is  the  largest  when  the 
molecule  is  next  to  the  ion  and  becomes  smaller  when  it  is  in  the  next  two  shells.  The  values 
of  the  dipoles  on  the  water  molecules  is  larger  when  they  are  around  F',  again  due  to  the 
large  field  emanating  from  the  ion. 

To  see  how  the  hydrogen  bonding  pattern  changes  when  water  is  around  the  ion  in  the 
cluster  we  calculate  the  percentage  of  0,1,2,3,4  and  5  hydrogen  bonding  water  molecules  in 
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the  first  and  other  shells.  The  histograms  representing  these  percentages  are  shown  in 
Figures  7  and  8.  As  we  can  see  from  these  figures,  the  first  shell  water  molecules  are  non¬ 
hydrogen  bonded  in  the  F’fHjOn  clusters  with  n  up  to  five  waters;  the  one-hydrogen  bonded 
waters  are  dominant  in  the  clusters  v.dth  n  greater  than  6  but  less  than  14,  the  two-hydrogen 
bonded  water  molecules  is  the  dominant  mode  in  clusters  with  14  and  15  water  molecules. 
For  the  Cl  fHjOln  clusters  the  first  shell  water  molecules  are  not  hydrogen  bonded  only  in  the 
clusters  with  n  up  to  two  waters;  the  one-hydrogen  bonded  waters  are  dominant  in  the 
clusters  with  n  from  3  to  10;  the  one-hydrogen  bonded  molecules  are  found  as  often  as  two 
in  clusters  with  n=  11  and  12  waters;  finally  the  two-hydrogen  bonded  water  molecules  is  the 
dominant  mode  in  clusters  with  13  to  15  water  molecules.  For  the  second  shell  water 
molecules  the  difference  in  the  hydrogen  bonding  pattern  is  not  that  pronounced:  for  the  F' 
(HjOn  and  Cr(H20)„  clusters  with  n=2,3,...9  the  most  often  found  are  water  molecules  with 
one  hydrogen  bond  per  molecule.  Two  hydrogen  bonds  is  the  most  dominant  mode  for 
molecules  in  larger  size  clusters  with  F'  ion.  In  clusters  with  Cl  ion  two  hydrogen  bonding 
molecules  dominate  in  the  secondary  shells  when  n=  10,11,12;  for  clusters  with  n  =  13,14  and 
15  the  amount  of  water  molecules  with  two  hydrogen  bonds  per  water  is  the  same  as  with 
three  hydrogen  bonds.  This  again  confirms  that  the  hydrogen  bonding  network  in  water  is 
stronger  in  clusters  with  Cl'  compared  to  the  clusters  where  the  F'  ion  is  solvated. 

As  we  have  mentioned  it  in  Introduction,  the  electrostatic  stabilization  energies 
obtained  from  the  photodetachment  spectra  of  Br'(H20)„  and  I'dljOl^  (n  =  1,2...15)  clusters  can 
be  scaled  so  that  the  data  fall  on  a  common  curve.  That  means  that  Br'  and  I'  ions  are  in 
similar  environments  in  clusters  with  the  same  number  of  water  molecules.  From  our 
previous  discussions  on  Cl'  and  F'  we  have  to  conclude  that  since  the  F'  ion  is  solvated  in 
water  clusters  and  Cl'  ion  is  not,  the  stabilization  energies  of  these  two  ions  in  water  clusters 
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can  not  be  brought  to  the  same  curve  by  a  simple  scaling.  Indeed,  this  is  confirmed  in  Fig. 
9  where  we  display  the  calculated  values  of  the  electrostatic  stabilization  energies.  The 
detailed  description  of  how  such  calculations  are  done  is  given  in  our  previous  paper  [27].  Let 
us  note  here  that  the  values  of  are  experimentally  measvu’ed  for  Cl'dljOln  (n=l,2...7) 
clusters  (our  calculations  are  in  good  agreement  with  these  measurements),  while  no  data  are 
available  for  the  clusters  with  F'  ion.  Therefore  only  future  experimental  work  on 
photodetachment  spectra  for  larger  sized  clusters  with  Cl'  ion  and  on  clusters  with  F‘  ions  in 
water  will  show  how  accurate  are  our  calculations. 

Very  often  the  cluster  structmes  are  discussed  in  terms  of  minimum  energy  structures 
and  transitions  between  these  minimum  energy  structures.  For  small  atomic  clusters  this 
method  can  give  a  complete  picture  of  the  structure  and  d5mamics  in  clusters,  since  one  can 
find  all  local  minima  of  the  clusters  [34].  For  molecular  clusters  a  comprehensive  search  for 
all  local  minima  is  quite  impossible.  Nevertheless  it  is  very  interesting  to  know  how  the 
structure  of  the  cluster  in  its  global  energy  minimum  looks  like,  since  this  is  the  structure 
the  cluster  would  have  at  T = 0  K.  Comparison  of  this  structure  with  the  structures  at  T  ~250 
K  can  give  us  an  idea  how  important  the  energy  and  the  entropy  components  of  the  free 
energy  are  in  determining  the  cluster  structures  at  T~250  K.  To  find  the  local  energy 
minima  for  some  of  OUT  clusters  we  quenched  500  structures  from  molecular  dynamics 
trajectories  of  the  clusters  C1(H20)„  and  F'(H20)n  (n= 6,7,8).  The  starting  points  for 
quenching  were  separated  by  2  ps  each.  The  quenching  was  performed  using  the  conjugated 
gradient  method.  In  Figure  10  we  represent  the  lowest  energy  structures  that  we  obtained 
from  our  quenching  procedures.  To  understand  the  geometry  of  these  clusters,  we  have  to 
remember  that  the  minimum  energy  structure  of  a  cluster  with  eight  water  molecules  is  a 
cube  [35].  As  we  can  see  from  Figure  10,  the  structure  of  a  cluster  C1(H20)7  is  also  a  cube, 
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where  Cl'  substituted  one  of  the  water  molecules  at  the  vertex  of  the  cube.  We  can  see  also 
that  the  structure  of  other  clusters  can  be  also  understood  as  being  derivatives  of  the  basic 
cubic  structure.  As  always,  it  is  hard  to  know  if  our  structures  represent  the  global  minima, 
but  the  similarity  of  the  structures  to  the  cubic  structure  of  the  established  global  minimum 
indicate  that  om*  structures  are  close  to  the  global  minimum  structures.  Fig.  10  also  shows 
that  both  P  and  Cl"  ions  prefer  to  be  on  the  surface  of  the  minimum  energy  clusters  with  6, 
7  and  8  waters.  But  as  our  molecular  dynamics  simulations  show,  at  T~250  K  the  F'  ion  is 
fovmd  inside  the  clusters,  while  Cl'  still  remains  on  the  surface  (see  Fig.  5a  and  5b).  That 
means  that  at  least  in  clusters  with  6,7  and  8  waters  the  solvation  of  the  F'  ion  is  due  to  the 
entropy  effect.  A  more  detailed  study  of  minimum  energy  structures  of  water  clusters  with 
anions  is  now  underway. 

CONCLUSIONS. 

We  have  performed  molecular  dynamics  calculations  on  Cr(H20)n  and  F'lHgOln 
(n=2,3,...15)  clusters.  These  calculations  show  that  the  F'  ion  is  solvated  in  these  clusters, 
while  Cr  remains  attached  to  the  water  in  the  clusters.  In  our  previous  work  we  used  the 
argument  of  similarity  to  show  that  the  I'  and  Br'  ions  are  also  attached  to  the  water  in  the 
clusters.  Very  recently  a  direct  molecular  dynamics  calculation  off  in  water  clusters  showed 
that  the  I'  ion  is  on  the  surface  of  the  cluster  [36].  Our  calculations  of  the  clusters  with  Br 
ion  also  show  that  the  ion  is  on  the  surface  of  the  cluster  [37]. 

Based  on  our  calculations  presented  above  and  on  our  previous  calculations  we 
conclude  that  F'  ion  is  solvated  in  clusters  with  up  to  15  water  molecules,  while  larger  size 
halide  anions  are  not  solvated  in  this  size  clusters.  Why  then  is  F'  solvated  in  water  clusters 
while  larger  size  halide  anions  are  just  attached?  To  answer  this  question  let  us  assume  that 
the  water  molecules  in  the  cluster  do  not  want  to  disrupt  their  hydrogen  bonded  network. 
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One  needs  a  rather  strong  perturbation  to  do  this.  Such  perturbation  may  be  produced  by 
the  F"  ion,  because  it  is  smaller  in  size  than  the  other  halogen  ions  and  therefore  the  electric 
field  around  it  has  a  larger  value.  But  even  then  the  perturbation  may  not  be  strong  enough 
(as  we  saw  from  minimum  energy  structures  for  PdijOln  (n= 6,7,8))  and  therefore  the 
solvation  of  the  ion  is  due  to  the  entropy  effect. 

One  can  ask  a  question  related  to  the  reliabUity  of  oiu*  potential  functions  used  in  the 
simulations.  As  we  have  demonstrated  above  the  calculations  with  these  functions  reproduce 
rather  well  enthalpies  of  formation  of  small  F  (H20)n  and  C1(H20)„  clusters.  Moreover  the 
electrostatic  stabilization  energies  in  small  Cr(H20)„  (n=  1,2,. ..7)  clusters  are  also  reproduced. 
The  electrostatic  stabilization  energies  of  Br  and  V  ions  that  we  can  calculate  from  the 
molecular  dynamics  simulations  can  be  compared  with  the  experiment  for  clusters  with  up 
to  15  water  molecules.  Unfortunately  this  is  not  the  case  for  the  clusters  with  C1‘  and  F'  ions. 
Therefore,  more  experimental  data  are  needed  to  confirm  our  conclusions  about  the  behavior 
of  Cr  and  F‘  ions  in  water  clusters. 
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Table.  II.  Calculated  (MD)  and  experimental  (exp)  heats  of  formation  of  X'CHjOln  clusters  in 
kcal/mol. 


n 

Cr  (MD) 

cr  (exp)^ 

F-  (MD) 

F'  (exp)^ 
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1.  Values  are  the  best  estimates. 

2.  Reference  [4]. 
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Figure  Legends: 

Fig.  la.  The  distribution  of  distances  from  the  ion  to  the  center  of  mass  of  the  cluster .  Solid 

line  is  for  the  CITHjOIt  cluster;  dotted  for  the  FCHjO)^  cluster 

Fig.  lb.  The  same  as  la,  only  for  the  clusters  with  15  water  molecules. 

Fig.  2a.  The  ion-water  oxygen  radial  distribution  fimction  for  F(H20)„  clusters.  Dashed  line 
is  for  the  cluster  with  n=3,  dotted  line  is  for  n=7,  solid  line  is  for  n=  11  and  dot-dashed  line 
is  for  n=15. 

Fig.  2b.  The  same  as  Fig.2a,  oiJy  for  Cr(H20)„  clusters. 

Fig.  3a.  The  running  coordination  number  for  the  F  (H20)„  (n  =  2,3,...15)  clusters. 

Fig.  3b.  The  same  as  Fig.  3a,  only  for  ClXHjOln  clusters. 

Fig.  4a.  The  angular  distribution  function  for  F  (H20)„  clusters.  Dashed  line  is  for  the  cluster 
with  n=3,  dotted  line  is  for  n  =  7,  solid  line  is  for  n=  11  and  dot-dashed  line  is  for  n=  15. 

Fig.  4b.  The  same  as  Fig.  4a,  only  for  Cl  (HjOn  clusters. 

Fig.  5  A  snapshot  from  a  tr^ectoiy  for  the  F'lHjO);  cluster  (top)  and  Cl  CHjO);  cluster 
(bottom). 

Fig.  6a.  The  values  of  the  total  dipole  moment  of  water  molecules  for  the  waters  in  the  first 
solvation  shell  (circles)  second  solvation  shell  (squares)  and  third  solvation  shell  (triangles) 
in  the  F  (H20)n  (n=2,3,...15)  clusters. 

Fig.  6b.  The  same  as  Fig.  6a,  only  for  CKHjOln  clusters. 

Fig.  7a.  The  percentage  of  0,  1,  2,  3,  4,  5  hydrogen  bonded  water  molecules  in  F  (H20)„ 
clusters.  The  water  molecules  belong  to  the  first  hydration  shell  around  the  ion.  We  used 
the  geometric  criteria  (see  reference  [38])  to  determine  the  presence  of  the  hydrogen 
bonding. 

Fig.  7b.  The  same  as  Fig.  7a,  only  for  Cr(H20)n  clusters. 
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Fig.  8a.  The  percentage  of  0,  1,  2,  3,  4,  5  hydrogen  bonded  water  molecules  in  F  (H20)„ 
clusters.  The  water  molecules  belong  to  the  second  and  third  hydration  shells  around  the  ion 
Fig.  8b.  The  same  as  Fig.  8a,  only  for  Cl  (H20)„  clusters. 

Fig.  9.  Electrostatic  stabUization  energies  for  the  Cl  (H20)„and  F  (H20)n  (n  =  2,3,...15)  clusters. 
Fig.  10.  F  (H20)„  and  CKHjOln  (n  =  6,7,8)  clusters  in  their  minimum  energy  configurations  . 
On  the  left  are  the  clusters  with  the  F'  ion  and  on  the  right  are  clusters  with  the  Cl  ion.  The 
upper  row  depicts  clusters  with  n  =  6,  the  middle  row  n=7,  the  lower  row  n  =  8. 
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